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SUMMARY 

A rigorous treatment is given of the problem of deter- 
mining the two-dimensional potential flow around 
arbitrary biplane cellules. The analysis invoices the 
use of elliptic functions and is sufficiently general to 
■include the effects of such elements as the section shapes, 
the chord ratio, gap, stagger, and decalage, which elements 
may be specified arbitrarily. The flow problem is re- 
solved by making use of the methods of conformal repre- 
sentation. Thus the solution of the problem of trans- 
forming conformaUy two arbitrary contours into two 
circles is expressed by a pair of simultaneous integral 
equations, for which a method of numerical solution is 
outlined. It is pointed out that 'an inverse method of 
transforming conformaUy two circles into the wing pro- 
files of a biplane arrangement leads readily to the devel- 
opment of relatedfamilies of biplane combinations. Flow 
formulas are developed giving the velocity and pressure 
at any point of the surface of either profile of the arbi- 
trary biplane arrangement, for any angle of attack. The 
theory of the monoplane wing section in potential flow is 
shown to be a degenerate case in which the elliptic func- 
tions reduce to trigonometric functions. The general 
method presented may be employed to determine the 
potential flow in any doubly connected region and hence 
may be applied to the single slotted wing or to the auxiliary- 
airfoil wing. 

As an example of the numerical process, the pressure 
distribution over certain arrangements of the N. A. G. A. 
4412 airfoil in biplane combinations is presented and 
compared with the monoplane pressure distribution. 

INTRODUCTION 

It is the purpose of this paper to develop a general 
theory of arbitrary biplane cellules of infinite span in 
potential flow. No attempt is made here to treat the 
case of finite span or to consider viscosity; rather it is 
the object of this work to bring the two-dimensional 
theory of biplane cellules in uniform, steady potential 
flow to the same degree of exactness and generality to 
which the two-dimensional monoplane airfoil theory 
has been brought. The analysis will be sufficiently 
general to include such elements as profile shapes, 
chord ratio, gap/chord, stagger, and decalage, and 



will contain as special cases the monoplane theory, as 
well as the theories of the slotted monoplane wing, of 
the auxfliary-airfoil wing, and of the influence of the 
ground or plane barriers on a monoplane airfoil in 
two-dimensional potential flow. 

In order to arrive in a natural manner at a perspec- 
tive of the biplane analysis it is advantageous to con- 
sider first the simpler case of the monoplane wing 
section and to keep in view the essential concepts that 
carry over to the biplane analysis. It is well known 
that by virtue of the methods of conformal represen- 
tation the two-dimensional potential flow around a 
single obstacle can be obtained by the following 
process. In the first place, a standard contour is 
selected, the region about which is simply connected 
and the flow function of which in uniform potential 
flow is known or obtainable. The transformation 
must then be found that transforms conformaUy the 
region of the given obstacle into this standard region. 
This transformation, in combination with the known 
flow function, gives the desired flow function for the 
obstacle. In the case of monoplane wing profiles, the 
standard flow region may be chosen to be that about 
a circle and the theorem which states that it is possible 
to transform conformaUy the contour of the given 
obstacle into a circle is known as Kiemann's theorem. 
(Of., for example, reference 1.) In the case of two 
obstacles, the region is termed "doubly connected" 
and the process is again applicable except that the 
standard doubly connected region is chosen to be the 
region about two circles. The theorem that states 
the existence of a transformation function bringing 
the doubly connected region (region of the biplane 
contours) into the region of two circles is Koebe's 
theorem (reference 2). 

The flow function giving the uniform potential 
flow for a circular cylinder is weU known and, in 
determining the flow about a monoplane airfofl section, 
the main problem is the transforming of the airfoil 
contour into a circle. In order to attain this result 
in a simple manner it is necessary to perform a few 
intermediate transformations. The airfoil profile it- 
self may be regarded as a contour described about a 
conveniently chosen line segment or chord. An 
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initial transformation of a simple type exists (the 
so-called "Joukowsky transformation") that trans- 
forms the chord into a unit circle and automatically 
maps the airfoil contour into a nearly circular contour 
described about the unit circle. There remains then 
the final task of transforming the nearly circular 
contour into a true circle, and this may be performed 
by a method given by Theodorsen (reference 3). This 
method leads directly to a simple integral equation 
which can be solved by a process of iteration or suc- 
cessive approximations and which converges with 
extreme rapidity. (Cf. reference 1.) It is important 
in regard to practical considerations to observe that 
the method is so powerful that one step in the process 
is quite sufficient in all ordinary cases. 

The standard doubly connected region has been 
chosen as the region about two circles; and it is worthy 
of mention that only as recently as 1929 was the 
complex flow potential for two circles rigorously 
developed (by Lagally, reference 4). Dupont, Bonder, 
and Muller (references 5, 6, and 7) have also con- 
tributed to this problem but Lagally 's solution is 
the more elegant. The flow function for two circles 
being known, the main problem in finding the flow 
about a biplane arrangement is the obtaining of the 
transformation mapping the two contours into two 
circles. 

In a manner analogous to the case of the single air- 
foil section, the contours of a biplane arrangement may 
be considered to be described about a skeleton of two 
conveniently placed mean lines or chords. Hence, to 
maintain the analogy it is seen that initially it is desired 
to find a transformation function which transforms 
the two line segments into two circles. This prob- 
lem has been touched upon by Kutta who has given the 
uniform potential flow for the special case of two paral- 
lel equal line segments (reference 8). The transforma- 
tion function bringing two circles into any two ■parallel 
line segments has been developed by C. Ferrari (ref- 
erence 9). In the first part of the present paper the 
more general problem of decalage of the line segments 
has been studied and a function developed that trans- 
forms two circles into any two nonintersecting line seg- 
ments in any relative positions. This function that 
transforms the skeleton or chords of the biplane ar- 
rangement into two circles also transforms the contours 
themselves into two nearly circular contours described 
about the skeleton circles. There remains then the 
problem of transforming the two nearly circular con- 
tours into two true circles. In order to accomplish 
this task, the method of Theodorsen is generalized in 
the present paper to apply to doubly connected regions 
by employing the concentric circular ring region as a 
standard region and by utilizing a Laurent series devel- 
opment instead of a one-way power series. There is 
obtained finally a pair of simultaneous integral equa- 



tions expressing the conformal representation of the 
two nearly circular contours into two circles. Just 
as in the case of the single integral equation in the 
monoplane case, there exists an analogous process of 
successive approximations or iteration that converges 
with the same remarkable degree of rapidity. 

The general transformation from the biplane con- 
tours to two circles together with the Lagally formula 
for the flow about two circles yields an expression for 
the velocity and pressure at each point of the surface of 
either profile of the biplane arrangement. There are 
two arbitrary circulations in the flow formula, viz, the 
separate circulations around each contour, and these 
are determined uniquely by applying the well-known 
Kutta-Joukowsky condition to the trailing edges of both 
contours, specifying thereby that the flow leaves these 
edges smoothly. 

In the case of monoplane wing theory it has been 
shown (reference 1) that theoretical shapes can be 
conveniently developed by an inverse method of 
" transforming conformally a circle into a wing profile. 
The Joukowsky airfoils and the other so-called "theo- 
retical" airfoils are special examples of this process. 
In an analogous manner it is possible to develop theo- 
retical biplane combinations by an inverse process of 
transforming two circles into two contours resembling 
wing profiles. A general and flexible method of obtain- 
ing these shapes is presented; the results are especially 
instructive in that, in this process, the integral equa- 
tions referred to in a preceding paragraph reduce to 
definite integrals. 

Elliptic functions arise in a natural manner in the 
analysis and the problem treated provides a good 
illustration of the power and beauty of these remark- 
able functions. The general theory of the single, 
arbitrary wing section is shown to be a degenerate 
case in which the imaginary period of the doubly 
periodic elliptic functions becomes infinite, and hence 
the elliptic functions reduce to ordinary trigonometric 
functions. A few pages are devoted to the monoplane 
theory in view of the light that it throws on the more 
general biplane analysis. 

Numerical results are presented only to furnish an 
illustration of the theory. In particular, the pressure 
distribution is determined for certain arrangements of 
the N. A. C. A. 4412 airfoil in biplane combinations. 
The elliptic functions that arise in the analysis and 
that are to be evaluated in a numerical case may 
fortunately, when necessary, be developed in rapidly 
convergent expansions. 

Statement of the problem. — The problem treated in 
this paper may be restated as follows. Given, an 
arbitrary biplane arrangement oriented in a specified 
manner in a nonviscous, incompressible fluid medium 
and translated with uniform velocity V. To determine 
the velocity and pressure distribution in two-dim en- 
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sional potential flow in the field of motion for all 
angles of attack, particularly, at each point of the 
surface of the biplane profiles. 

As has been pointed out, it is well recognized that the 
aforementioned problem may be treated in two stages. 
In the first place, the complex function expressing the 
conformal transformation of the region of the biplane 
into a standard doubly connected region must be 
obtained and, finally, the complex flow function for 
this standard region, which is chosen as the region 
about two circles, must be known. The region ex- 
ternal to the two contours of a biplane arrangement 
will be brought into the region about two circles by the 
intermediate use of two nearly circular contours.- 
Before this result can be accomplished, however, it is 
desirable to discuss several preliminary transforma- 
tions. 

1. PRELIMINARY TRANSFORMATIONS 

First, the transformation bringing the region external 
to two nonintersecting circles (t plane) into the annular 
region between two concentric circles (w plane) will be 
obtained. This annular region will then be mapped 
into a rectangular region (s plane) and the rectangular 
region into the region about two line segments (u 
plane). (See fig. 1.) 



t plane 



The curves r 2 /j*i= constant are circles with centers 
lying along the t 3 axis (theorem of Apollonius). This 
family of circles contains the points Qi and Q 3 as 
limiting circles of zero radius. For points in the 
upper half plane (i 2 >0) we have r 2 /r!>l, for points in 
the lower half plane r 2 /rj<l, while on the t t axis, 
rifr i=l. The curves given by y 3 — -y 1 =constant also 
form a family of circles (theorem of the constant angle 
subtended by the chord of a circle) which is orthogonal 
to the first system and each circle of which contains 
the limit points Q x and Q 3 on its circumference. 

A new complex variable w—ce 1 ^ 13 is now introduced 
by the following relation 



Hence 



Also 



C t—lC Ti 



M=log ri 



t=ic 



w+c 



w—c 



(2) 



(3) 



(4) 





These equations transform conform ally the coaxial 
system of circles of the t plane into a concentric system 
of circles in the w plane. In particular, two circles Kj 
and K s in the t plane, Ki located in the upper half 



s plane 



1— A 




(a) (b) (c) (d) 

Fioube 1.— Mapping oh (a) two coaxial circles In the t plane Into (b) two concentrio circles In the w plane, (o) rectangular region In the s plane, (d) two line 

segments in the u plane. 



Transformation of a coaxial system of circles into a 
conoentrio system. — A coaxial system of circles may 
be described most simply by the use of bipolar coordi- 
nates. Consider a complex t plane where f=fi+#s- 
Let Qi (0, ic) and Q 3 (0,—ic) located on the U axis be 
the origins of two polar coordinate systems r u yi and 
f2>73- The variable t may be written in the two forms: 

f =ic+rie n 'i= — ic-\-r s e"i 
Then in the relation 



t+ic r a 



(1) 



t—ic r x 

there are expressed in a convenient form, the bipolar 
coordinates p and -y 2 — yu (See fig. 1 (a).) 



plane and K s in the lower half plane and defined by 
logr s /ri=aandlogr 2 /ri=— 0, respectively (a>0, /3>0), 
transform into two concentric circles B x and B 2 about 
the origin in the w plane, of radii ce" and ce - ", respec- 
tively (fig. 1 (b)). It is noted also that the h axis 
transforms into the circle of radius c in the w plane, 
and that the region at infinity in the t plane maps into 
the neighborhood of the point w=c. It may also be 
remarked that the circles orthogonal to Ki and K 2 
transform into radial lines through the origin. 

Transformation of the circular systems into rectan- 
gular systems. — There is now introduced another vari- 
able «=X+w defined by the relation 



. , to . , t+ic 



(5) 
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Separating into real and imaginary parts 

X=— (72— Ti)=— 0 

1 r a (6) 
v=log-= M 

Hence the variable s may hereafter be denoted by 

s= — d+i/s. or also by s=X+ijtt 
Also from (5) we have 



— - rp~ is 



w=ce' 



and 



(7) 
(8) 



The circles in the t plane, r 2 /r!= constant (or the 
circles in the w plane ce' 1 =constant), correspond 
uniquely to the straight lines p= constant in the s 
plane. In particular, the limiting points Qi and Q% 
correspond to /i=a> and /i= — °o, respectively. Also 
the t t axis corresponds to the axis /x=0, the point at 
infinity in the t plane going into the origin 8=0. 
The circular arcs 7[— 72= constant between Qi and Q 3 
correspond to the lines X=constant. It is noted, how- 
ever, that this latter correspondence is infinitely many- 
valued since the addition of integral multiples of 2x to 
7i or 7a does not alter the circular arc considered; 
hence, 71— 7 2 =constant corresponds to the infinite 
number of parallel lines X==constant+2ta-, where k 
is any integer (fig. 1 (c)). 

The whole t plane has thus infinitely many values 
on the s plane, but there is a one-to-one correspondence 
between the whole t plane and a strip of width 2v 
bounded by two parallels to the p axis. In the fol- 
lowing investigation, the strip in the s plane bounded 
by the lines X=— 7r and X=x will be considered as 
the representation of the t plane cut along the length 

Qi Q*. 

Equation (5) thus defines a conform al transforma- 
tion of the coaxial system of circles in the t plane, or 
of the concentric system of circles in the w plane, to a 
rectangular system in the s plane. In particular, con- 
sider again the two" definite circles Ki and K 3 of the 
coaxial pencil. The circle K t is defined by log r 2 /r a = 
11= a and the circle K 3 by log r 2 /ri=n=—fi where a and 
P are positive constants. It is then noted that the 
region of the t plane external to the circles Ki and K 3 
(or the ring region within B 1 and B 2 ) corresponds 
uniquely to the rectangular region bounded by the 
lines n=a, n= — /3, X=— -k, and \=7r. (The two sides 
X=— tt and X=7r correspond to the right and left 
edges respectively of a cut along the t 3 axis drawn 
between the two circles.) The rectangle contains, 
necessarily, the point s=0 as an internal point. 

Geometrical relations. — Attention may be momen- 
tarily diverted to some geometrical relations existing 
in the various planes. Let the radii of Ki and K 2 be 
a and 6, respectively, and let the centers of Ky and K 3 
be situated at O t and 0 2 , respectively (fig. 1 (a)). 



The quantities a and b may be expressed in terms of 
a and /S. The equation of Ki in bipolar coordinates 
is, by equation (2) 



ri~ 



t+ic 
t—ic 



=e a 



Writing t=ti+it 2 there results upon expansion 

h s +t 3 s —2ct s coth a+<?=0 

which is the equation of a circle whose center Oi is 
situated at 

fa=c coth a 

and whose radius is 

a=c csch a 

Similarly, for the second circle K 3 , the center 0 2 is at 
tt=—c coth 0 

and the radius is 

b=c csch & 

Denoting by d the center-to-center distance O1O3 
(fig. 1), there may be written the equations: 

a=c csch a 
b=c csch j3 
d=c (coth a+coth /3) 

which suffice to fix a, b, and d in terms of a, /S, and c. 
Forming the auxiliary quantity (P—aP—b*, it is found 
thattP— a*— 6*=2a5cosh(o+/3) It immediately follows 
that 

c=-j sinh (a+/5) 
b 



(9) 



sinh a=^ sinh (a+/S) 



sinh /3=2 sinh (a+/3) 



(10) 



We observe that the quantity a-t-jS on the right-hand 
side is expressed in terms of a, b, and d by the relation 

cP-a'-b* 
' 2ab~ 



cosh (a+/3) = 



2. TRANSFORMATION OF TWO CIRCLES INTO TWO 
ARBITRARY LINE SEGMENTS 

The transformation that maps the rectangular region 
in the s plane into the region external to two noninter- 
secting line segments in a u plane will now be derived. 
In combination with the preliminary transformations 
of the preceding section this result will then transform 
the region of the two circles Ki and K 3 of the t plane, 
or the ring region of the w plane, into the region of 
the line segments. Let the u plane (fig. 1 (d)) contain 
the two line segments c 2 and c a and let u(t)=X+iY be 
the analytic function that transforms the circles Ki 
and K 2 into the desired line segments. With no loss 
in generality, the system of coordinates in the 11 plane 
may be so chosen that the X axis is parallel to c 2 . Let 
the line segment c x be inclined at an angle —5/2 with 
respect to the X axis. (The negative sign before 5 
is a matter of later convenience; fig. 1(d)) may be 
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regarded as illustrating the definition of positive 
decalage.) Let h and Z 2 denote the two lines fi=a 
and n=—p of the rectangular s plane that corre- 
spond to Ki and K 2 , then, it is evident that 

dY S 

for points of K 2 or Z 3 and jy=— tan ^ for points along 

Ki or li. Let/(s)=^ be the derivative of the func- 
tion u(s) that gives the desired correspondence between 
the u and s planes. From the well-known property of 
conformal mapping, viz, that tangents at correspond- 
ing points in the two planes differ in direction by the 
argument of the derivative function, it follows that 
the argument of f(s) equals 0 (or n) along Z 2 and equals 

—8/2 (or — 5/2+H-) along l u or 

/(s) is a real quantity along l 3 
f(s)e K/2 is a real quantity along li 

By a principle of Schwarz the as yet undetermined 
function /(s) has the property of being extended by 
analytic continuation to the whole strip region in the 
s plane (fig. 1(c)) for, since /(«) is real along Z s , its 
values for a pair of reflected points mirrored in the 
line Z a are conjugate complex. Similarly the function 
/(s)*"' 3 may be reflected about the line h. With 
successive alternate reflections in Zi and Z a /(«) takes 
on values as shown in figure 2. For every two succes- 
sive reflections the original values of /(«) are repeated 
except for a multiplying factor e~ a . Hence it is clear 
that /(s) must satisfy the relation 



(1) 



Also, since /(s) is a single-valued function of t, it 
satisfies the condition 



%+2^)=/( 5 ) 



(2) 



If 5=0, then e~ is —l and it is seen at once that/(s) is 
then a doubly periodic function, hence an elliptic func- 
tion (of the first kind), of real period 2w=2ir and of 
imaginary period 2<a'=2i(a+p). In the general case 
where 3p*0, the function /(s) is not a purely doubly 
periodic function but, since one of its periods gives 
rise to a multiplying factor, is an elliptic function of 
the second kind. It is completely determined, except 
for a constant, by its behavior at its poles, in the 
neighborhood of which the function becomes infinite 
(Hermite's theorem). In the present analysis, we 
shall consider the fundamental periodic rectangle as 
formed by the original transformed rectangle and its 
reflection in the line Z 2 (fig. 2). 

We now investigate the poles of the function 
dtL 

/(s)=^-- We assume that at infinity MH4 ^ order 
that the regions at infinity in the u and t planes be 



equally magnified and map into each other (except 
for a possible change in direction), and we have 



\du\ 

|2fW 



= 1 



or noting equation (1.8) 1 
du 



du _ ds I 2c _j .g. 



■flsje 



v , v//y. 



3s 



2jr- 



s plane 

— m 

T A 

_- -jW 

.—-/We" 

Figure 2.— Illustrating the properties of ^jr/M In the strip region: AO Is the 
conjugate complex quantity of Ks). 

This relation shows at once that/(s) possesses a singu- 
lar point at 8=0, and hence has one also at the point 
obtained by reflection of s=0 in Z 2 , viz, s= — 2ij3. In 
the neighborhood of these points /(s) becomes infinite 

in the order of f as f-»oo or as ^ as s-»0, i. e., has 

poles of order two at the origin s=0 and at s=— 2i/?. 

i This notation denotes equation (8) of sec 1. 
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The fundamental function having a single pole of 
first order (with residue unity) at the origin and satis- 
fying the foregoing period requirements is (reference 
10, p. 416, and reference 11, p. 369) 



A(s)=- 



-ts 

' <r{8+S) 



(4) 



where a denotes the sigma function of Wei era trass and 
possesses the following period properties 

o-(u+2co) — e-*' ( "+«><r(i0 
o-(u+2w') = — e-V(«+-o 0 -( u ) 

The expression (4) for A(s) may also be written as 
follows 

H'(Q)H(8+8) 

- ~nmm~ 



(5) 



(4') 



where the Jacobi IT(eta) function is defined by the 
equation (cf. references 11 and 12) 



ff(u)=2 2 '/ 1 sing-22»/*8in^+2^sin^- 

and possesses the period properties 
H{v+2u) =-E(u) 



(6) 



where 



fl"(ti+2«')=-2- 1 e ' Z>1 E(u) 



q_=e 



The relation existing between the <r and H functions is 
(reference 11, p. 488) 



<r(u)=e ' 



i7U» 



(7) 



A function such as we are seeking, having a single 
pole of the second order at the origin and the required 
period properties, may be obtained by taking the 
negative derivative of A(s) with regard to s. From 
equation (4') we have 



A'{8)- 



<M(s) 



ds 



ff'(0) d m 



The function /(s) is now determined except for con- 
stants ai and Oa is given by 



du 



/(s) =^=ai4' («) +a>A' (s+2ip)+a 3 



(9) 



To determine the constants, observe that by means of 
equation (3), and by the fact that the expansion of 

A' (s) about the origin begins with the term p> we have 
in the neighborhood of s=0 



|arl«_. 



B =] 



and since from equation (1.8) |f| t _» = 



there re- 



sults |ai|=2c. Thus the magnitude of ai is determined 



and, in general, we may put a i =2ce t '' where 7 is an 
arbitrary real parameter that determines the stagger 
of the segments, and the significance of which will be 
seen shortly. It may be observed at this point that 
with ffi 1 =2ce < T the following relation holds 



~dT\t. 



(10) 



i. e., the regions at infinity in the u and f planes agree in 
magnitude but differ by angle 7 in direction. In 
order to determine a 3 it is sufficient to recall that/(s) 
must remain real on U hence it may be seen that a? 
must equal 2ce~ n . Then finally equation (9) may be 
expressed as 

^=2c[A'(8)e t ->+A'(8+2ip)e- i - f ] (11) 

The general function relating the u and s planes is then 
by integration ■ ith regard to s 

v,(8)=-2c[A(s)e»+A(8+2ip)e- {y ]+k (12) 

where the function A(s) is given by (4) or (4') and 
where k is an arbitrary constant that is independent of 
s but may contain the parameter 5. s 

The singular points of transformation (12) are given 

by the roots of the equation ^=0. It is possible to 

draw at once certain conclusions with regard to the 
singular points. There exists a theorem on elliptic 
functions (reference 11, p. 366) which states that the 
number of zeros of an elliptic function (of the first 
or second kind) in a periodic rectangle is equal to the 

number of poles. Since /(s)=^ has 2 poles of second 

order in the periodic rectangle it follows that the 

equation ^=0 possesses 4 roots in this rectangle. It 

is demonstrable without difficulty that 2 of the zeros 
are located on the boundary n=a and the remaining 
2 on the boundary /x= — jS. These zeros correspond 
to the end points of the line segments Ci and c 2 of the 11 
plane. It may be stated for reference that on l 3 the 
singular values of X are obtained from the equation 



du 

d~8 



=0=A' (\—ip)e iy +A' (X-H/3)e- 



or since A'(X— ijS) and A'(\+ip) are conjugate com- 
plex quantities, the singular points are given by the 
solutions of 

Re. ^'(X-;/3)e ,T =0 

where Re. denotes "real part of." On l u we employ 
the period property (1) of A(s) and obtain for tho 
equation satisfied by the singular values of X 

Re. ^l , (X+ia)e <cr+8ZS5 =0 



> The remainder of this section is commentary to this equation. The reader 
may, without loss of continuity, proceed to seo. 3, p. 66. 
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Developments of A(s) convenient for numerical pur- 
poses will be discussed shortly. It may be of value to 
consider first several useful special examples of equa- 
tion (12). 

(a) Parallel segments of zero stagger 3 (5=0, 
7=0). — It is necessary to observe first the limiting 
form of A(s) as 5-K). From equations (4) or (4') (or 
cf. reference 10, p. 425) we have that 



s(^«-i)-7w-?-«"-; 



Iim 
s 



(13) 



where the various forms are equivalent. The func- 
tions <r and f are Weierstrass elliptic functions, and 
H and Zi (eta and zeta functions) are the elliptic 
functions of Jacobi and Hermite. Then putting for 

4c 

convenience the arbitrary constant k=- — 1 



■y— ic m equa- 



tion (12) we obtain 



t plane 



or we may put down the complete equation for refer- 
ence as follows. Noting that 

z.'w— IDw-^w+f 

where the Weierstrass p function is denned by ■ 
2»(s)=— ^f(s), and writing s=\+ia the equation 
determining the singular points of c x is 



IT 



(15) 



The addition theorem (reference 10, p. 140) of the p 
function may be written 

Here, p(a)=— p(ia) and the bar designates that the 
elliptic function p is based on periods 2u and 2o>' cod- 




c, 
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Fioubb 3.— Illustrating the case of parallel line segments: (a) t plane, (b) Ui plane ({-=0, 7=0), (0) Dj plane (J-=>0, t=i/0, (d) U plane (8=0, 1 arbitrary). 



(14) 



= -2c[Z 1 (s)+Z 1 (8+2ip)]—ic 
or (as given by Ferrari), 

= -2c[r («) +^s+m -^]+const. 

Observing that the Z y function has the following period 
properties 

Z 1 (8+2o>)—Z 1 (.s)=0 
Z 1 (s+2o>')-Z 1 (8)=-i 

it follows that for y.=a, y=c and for n= — /3, y= — c. 
Hence the gap Q between the line segments is 2c. 
The case is illustrated in figure 3 (b). 

The singular values of X are in this case given by 

Ee. Z/(s+ia)=0 
Ke. Zi'(s— i/3)=0 



The case of parallel line segments has been studied by Ferrari (reference 9). 



2 W ' 

jugate to the periods of p, i. e., 2&>=^- and 2&/=2i&> 

(reference 10, p. 32). Making use of (16) and of the 
differential equations for the p and p functions: 

p'W=^W-g s p(\)-g s 
p' (a) 2 =4i3 3 (a) -g 3 p (a) +g 3 

equation (15) becomes 

Ap(\y+Bp(\)+C=0 (17) 

where 

A=4[jj/7r— p(a)] 

B=4:p( a y+8(r,Mp_( a )-g s 

C=4( v /ir)p( a y+g s p( a )-£g 3 

This equation determines the two singular points 
for the upper segment c t . In general, there is only one 
positive value of the root p(\), hence the singular 
points, X=±X„ are symmetrical with respect to the 
origin. The negative root does not give real values 
for X. For the lower line segment c 2 it is only neces- 
sary to replace a by /?. (Cf. reference 10, p. 272: 
Given p(X), to find X.) 
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The line segments c x and Cj given by transformation 
14) are without "stagger" since the midpoint of each 
segment is located on the y axis (fig. 3 (b)). In order 
to obtain further insight into the general transforma- 
tion let us consider the case 5=0, 7=90°. 

(b) Tandem parallel segments ^5=0, Y=j0- — In 

tin's case let the arbitrary constant k=—c in trans- 
formation (12), and noting equation (13), we obtain 
Lim u(s)=Ui(s) 

3-0, 7-§ 

=—2ic[Z 1 (8)—Zi(a+2i0)]—c (18) 
=-2ic[r(s)-r(s+2i/3)]-c 

It can be shown directly that the singular values of X 
are 0 and ir, for we have 



and if 



p(u)=p(v) we must have 
u= ±»+2mu+2mV 



(19) 



where m and m' 



are any integers. Hence, writing s=X-Hm we have 
for the solutions of (19) in the fundamental periodic 
rectangle A,=0 and X,=x. Figure 3(c) shows a 
typical correspondence for this case. 

(c) Parallel segments of arbitrary stagger. — Let 

the arbitrary constant &=^r — ie^ cos y—c sin y in 

equation (12), and noting relation (13), 



Lim u(8)=U(j8)=Ui(j8) cos y+U 3 (s) sin y 
= — 2c(Z 1 (s)e"+Z 1 («)«-'»)— ice-" 



(20) 



The function U(s) is the general relation bringing the 
region about any two parallel line segments into a 
rectangular region in the s plane. In this transforma- 
tion the values of the parameters a, 0, and y suffice 
to fix uniquely the chord ratio Ci/<h, the gap/chord 
6/C3, and the stagger/chord S/c 3 of the parallel seg- 
ments. The gap between the line segments is 2c cos y. 
Figure 3(d) illustrates the definitions of the various 
quantities. In the general case of parallel segments 
of arbitrary stagger, there are given the three ratios 
Ci'.Cj'.G'.S and the parameters a, /3, and y are to be 
determined. This problem involves the solution of 
transcendental relations; in this connection it is con- 
venient to draw up charts, e. g., figure 4. This 
figure shows a cross plot that presents 6/c 3 , S/c t in 
terms of to' and y in the case of equal chords Ci/c 2 =l, 
i. e., a=/S. 

The singular points of equation (20) are defined by 
the relation 

^=0=p( 8 )en+p(8+2ip)e-"+^ cos y (21) 



Writing s=X+ia in (21) we have 

cos y\jp(\+ia)+p(\— ict)+2i}/ir] 
+i sin 7[p(X+ia)— p(\— ia)]=0 
Elmploying the notation of equation (17) and the rela- 
tions preceding equation (17), there results 

^ 4 (X)+a 3 p 3 (X)+a s ^(X)-l-a l3 >(X)-l-ao=0 (22) 




SO 40 60 

Angle of stagger, r .degrees 



80 



Fiottbs 4. — Ohart presenting the angle of stagger i and the ratio of the periods 
against gap/ohord 0/e», and stagger/chord Sla In the special case of equal 
parallel line segments, (S=0, a=P) 



where 



a i =A 3 <P 
a 3 =2ABd*—4 
a 1 =(B i +2AC)d 1 
a 1 =2BC(P+g 2 

a 0 =O 2 d 2 +g 3 and d= 



1 cot 7 



This equation suffices to determine the values of X 
corresponding to the end points of the upper line seg- 
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ments Ci. In the case of the lower segment a is to be 
replaced by /S. In general, there are only two positive 
solutions for p(\); and it may be observed that the 
solutions for both angles of stagger ±7 are contained 
in equation (22). (In sec. 4 approximations for the 
singular values of X are given by simple formulas.) 

Jacobi series — In order to obtain further insight 
into the general relation (12) and to separate u(s) = 
X-\-iY into its real and imaginary parts it is necessary 
to revert to Jacobi expansions. These developments 
are especially useful where numerical evaluations are 
required. 

By reference 10, page 416, we have the following 
expansion for the function A{a): 



A(8): 



-3s, 

e a <r(8+5) 
<r(«M«) 



=l(cot |+cot £)+2 S S a 5 """ sin (ms+nS) (23) 
where 

The expression for Z x («) occurring in the case of parallel 
segments is 

Um(A(8)-^ cot 0=Zi(s) 

=1 cot |+2 S S 2 2 *" 1 sin ma (23') 

To separate A(s) into real and imaginary parts, re- 
place s by X+itt and note that 



cot 



X-H/t _ sin X— i ainh y 
2 — cosh y— cos X 



(24) 



Then 



where 



A( S )=M(X, M )+#V(X,M)+^cotJ (25) 



M(X ^ )= ~2(coshT-cos X) 

CD CO 

+ S 2 a 2 ™" sin (mX+n5) cosh my 

pr f . . sinb y 

iV ^' ,t; ~ 2 (cosh yx-cos X) 

O* CO 

+ S S 2*"" cos (mX+w5) sinh my 

bj-1 n-1 



< Thovaldeofs-e-fat/n mayalv/aysbekeptle5sUianr^=0.0432byresortU]gwhen 
necessary, a. g., a+£ small, to transformations that interchange the real and Imagi- 
nary periods of the elliptio functions (reference 10, p. 266). Thus the expressions 
can always be made to converge very rapidly. Indeed, there exist several other 
expansions for A(fi) which thongh less simple in form are more rapidly convergent 
than the formula given here (reference 10, p. 422). 

136092—37 5 



Let us put the arbitrary constant k equal to 
(2c ,cot ^— ic)cos7— csin-y in equation (12) and sepa- 
rate u(s) as follows 
u(s)=X+iY 

=«i(s) cos y+ih(s) sin y (26) 

where 

u 1 {s)=x- s riy 



= -2c[A(s)+A(.s+2ip)]+2c cot ~ic (27) 



and 



= -2ic[A(s)-A(8+2ip)]-c (28) 
It is evident that Ui(s) and 1*2(5) are generalizations 
of Ui(s) and Z7 2 (s), given by equations (14) and (18) 
for the cases (5=0, 7=0) and (5=0, y=vj2), respec- 
tively. 

Employing relation (25), equation (27) giving u^s) 
may be separated into 

i=-2c[M(X, y)+M(\, y+2f})] 1 
y=-2c[N(\, y)+N(\, M+2/J)]-cJ {Z7&) 
It is observed that for 11= — /3, the coordinates become 

2te=— c J K J 

Equation (27a) is most useful in the neighborhood of 
H=— /3. Eor values of y near a, relation (27) is first 
rewritten by making use of the period property 
A(8+2o) , )=A(s)e~ u , and we have 

x=—2c[M(\, y)+M(\ y—2a) cos 5 
+N(\, it— 2a) sin 5] +c sin 5 
y=— 2c[iV(X, M)+iV(X, it — 2a) cos 5 
— M (X, tt— 2a) sin 5] +c cos 5 
For At=«, equation (27c) becomes 

X a =—2c[M(\ a)(l+COS S)—N{\, a) sin 5] 
+c sin 5 

2/a=— 2c[2V(X, a)(l-cos 5)— M(X, a) sin 5] 
+c cos 5 

It may be remarked that equations (27), which hold 
also for the special case 5=0, show immediately that 
in this case yp=—c and y a —c, or that the gap of the 
parallel segments is 2c. In the general case (5^0, 
7=0) it is clear from (27d) that the point (x a , y a ) 
= (0, c) lies on c x . • The "gap" as measured along the 
y axis (L. e., from x?=0 to x a =0) is therefore again 2c. 
The effect of decalage may be considered to a first 
order to be a rotation of the segment c t for the case 
(5=0, 7=0) by the angle 5/2 about the point (0, c). 

Employing relation (25), equation (28) giving u 3 (s) 
may be separated into 

x'=-2c[-N(\y)+N(\,y+2p)]-cl 
y'=-2c[M(\y)-M(\,y+2f})] ^ / ^ 8a; 
It is observed that for y=—fi the coordinates become 
i/=-4ciV(X,/3)-cl 

y/=o I 



(27c) 



(27d) 



(28b^ 
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(28c) 



(28d) 



In the neighborhood of /i=a, it may he preferable to 
express equation (28a) as follows 

x'=— 2c[-N(\ >f i.)+N(\,n— 2a) cos 

—M{\y.—2a) sin 5]+c cos 5 
y'= — 2c[M(\,n)— M{\,ix— 2a) cos 5 
— N(X,fi— 2a) sin 5]— c sin 5. 
For n=a the coordinates are seen to be 
x a '=2c[N(\a)(l+cos 5) 

+M(\a) sin 5]+c cos 5 
y a '=—2c[M{\,a) (1— cos 5) 

+N{\,a) sin 5]— c sin 5. 

For 5=0 it is clear from (28b) and (28d) that 
y a '=j//=0. In general, for 5^0 it can be seen that 
the coordinates (x</, y a ')=(—cfi) satisfy equation 
(28d), hence this point lies on c x . 

The general equation (12) or (26) may now be 
separated into 

u(8)=X+iY (29) 

where 

X=x cos y+x' sin y 

Y=y cos y+y' sin y 
In particular, it is clear from the foregoing that the 
lower segment c 2 is situated at Yp= —2c cos y and that 
the point {X a , Y a )=(—c sin y, c cos 7) lies on the 
segment Ci. 

If there are given any two line segments in position, 
the three ratios ^'.Cs'.OiS are known (in addition 
5 is known), and the quantities a, (), and 7 are to be 
determined. Equation (29) is transcendental and a 
direct solution for a given case is not available; 
however, an indirect procedure of building up charts 
similar to figure 4 (for which 5=0) for different values 
of 5 may be resorted to. The case of parallel segments, 
as well as the degenerate monoplane case (cf. sec. 
4), will prove helpful in this procedure. 

du 
3s 



For later reference, the derivative expression 
may be put down. We have 



(30) 



where 



du dui , din . 

=P(\v)+iQ(\n) 

±^=i[A'(8)-A'(s+S!if})] 

=P'(V)+*2'(V) 
In order to determine P, Q, P', and Q' , the following 
development is noted 



4 sin 1 s/2 m _, „, 
=M'M+iN'(.\n) 



S X! mtf™ cos (ms+n8) (31) 



where 

-v 1 — cos X cosh /1 
M (V)= 2(cosh m-cos X) 



OS 00 



-2E mg 2 " 1 cos (mX-fttS) cosh rati 



Tj/r\ \ — sin X sinh ft 

iV ^ ,/xj_ 2(cosh /i-cos X) 



OO CP 



— X) 2 ruf™ sin (mX+n5) sinh m/i 

m— 1 n— 1 

Hence, 

P=M'(X, M )+M'(X,^+2/3) 

Q=N'{\ri+N'(\p+2(}) 
P'=_j\P(X, JU )+i\P(X, /t +2/3) 

<2'=M'(X,,x)-M'(V+20) 
And finally 

die 

-^=2c[P cos 7+P' sin y+i(Q cos y+Q' stay)] 

(32) 

The equations of this section may be simplified in 
the noteworthy special case in which a=j8. The 
constant 2i/3 is in this case equal to half the imaginary 
period, i. e., 2ij8=w' and, in particular, the line seg- 
ments d and Cs are equal. By reference 10, page 422, 
we have 

=e _ ^| esc l+SSa**- 0 " sin[nw+(»-l/2) 5]) 

This expression may therefore replace A(s+2ifi) in 
equations (26), (27), and (28). Similarly in equation 
(30) A'(8+2ip) may be replaced by A'(s+w') where 

A'(8+u')= — e~% SS m2 l2n -« " cos [ms+ (»— 1/2)5] 

m=i n— 1 

3. TRANSFORMATION OF A NEARLY CIRCULAR RING 
REGION IN THE w PLANE INTO A TRULY CIRCULAR 
RING REGION IN THE z PLANE 

In the foregoing sections, there have been obtained 
the equations transforming the region external to two 
circles tin the t plane); or the annular region between 
two concentric circles tin the w plane) ; or also a rec- 
tangular region tin the s plane); into the region 
external to any two nonintersecting line segments (in 
the 11 plane). It may now be imagined, for definite- 
ness, that two airfoil profiles are generated about the 
two line segments as chords in the u plane (fig. 5 (a)). 
In the plane of the rectangle, the two profiles will 
correspond to curves of small amplitude extending 
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from X= — ir to X=7r near the boundary lines n—a 
and ii=—p, respectively. In the ring region, the 
profiles will correspond to two nearly circular contours 
forming an annular region (fig. 5 (c)). It is intended 
to show how this annular region may be transformed 
into a concentric circular ring region (fig. 5 (e)). 

At present it is assumed that the nearly circular 
ring region in the w plane corresponding to a given 
biplane cellule in the u plane is known. It is observed 
that this knowledge implies that equation (2.12) may 
be inverted and the variables X, y. solved for in terms 
of x and y. How this task may be done is taken up 
in section 4. It is recalled here that the variable X 
corresponds to — 0 (equation (1.6)) and that the 



where the radii are respectively, 

#1=0*% R 3 =ce?i, o-!>0, <r 2 <0 

(At times it will be found convenient to denote a x by 
a.' and tr 2 by — jS'). 

Let the function that transforms the w plane con- 
formally into the z plane be written as 

w=ze hW (5) 
where z=ce a+{ *=Re i * and where h(z) represents a 
Laurent series with complex coefficients: 

A(2)=a 0 +S(a»s n +a-„2-») (6) 
i 

where R 2 < \ z | <2?i, or with z=Re t * 

h{z)=f(R,<p)+ig(R ><P ) (6') 




u plane 




(a) 



(b) 









to plane 
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f m 
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z plane 












FiauBE B.— Mapping of: (a) two contours At and a of a biplane arrangement Into (b) two curved lines L% and Lt In the » plane, (c) two nearly circular contours 
Bi and flj of an annular region In the w plane, (d) two nearly circular contours JTi and Ki In the ( plane, (e) two true circles Ci and ft of the concentrlo ring region 
in the z plane. 



neighborhood of the point w=c, which corresponds to 
the region at infinity in the t or u planes, must be an 
internal point of the annular region. 

The annular region in the w plane then contains 
two boundary contours, an outer contour Bi and an 
inner contour B 2 . Let the contour B x be defined by 

w=C€Pt (0)+u (!) 
and the contour B 2 by 

w=ce»> un+u > (2) 

where the range of 0 may be chosen as O^0;£2n-. 

Consider now a z plane (fig. 5 (e)) containing two 
concentric circles about the origin, an outer circle d 
that corresponds to Bi and an inner circle 0 2 that 
corresponds to B 2 . The circle d may be defined by 



and the circle 0 2 by 



(3) 
(4) 



It is seen that on O u 
to 

log -=h(z)=h(R u <p) 

=f(Ri,<p)+ig(.Ri,<p)=Hi—<ri+i(d—<p) 1 
or, in short, 

=/i (<p) +Wi (<P) =m—<ri+i(0—<p)i (7) 
where (0 — <p)i means that the quantity 0 — <p is evalu- 
ated around d. On 0 2 similarly 

haiv) =f %(.¥>) +igi (<p)=ftt— ff 2 +i(d— <p) 3 
Let the complex coefficients in equation (6) 
expressed as 

a*=A n +iB n 
a- n =A-.„+iB-. 
Then, from equations (6) and (7), it is found that 

fib)=A 0 +ii[(AJZ l *+A^ l R 1 -') cos n? 

1 — (BjRi'—B-JRf*) sin n<p] (10) 
9i(<p)=B 0 +iO l [(BnRi n +B- x R 1 -*) cos n v 

+ (A n R 1 n -AJir*) win?] (11) 



J 



(8) 
be 

(9) 
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Similarly 

/ s ( SP )=^ 0 +2[(^nB 2 "+A_ H B 2 - B ) cos TUp 

— (BJ2 2 »— B_Bj-") minv] (12) 

ffa ( ¥ ,)=5 0 +S[(-B n B s B +.B_,5! s - B ) cos 

+ ( J 4„ J B 2 n -A. n B 2 - n ) sin Tip] (13) 

From equation (10), 

A_ B Br n =a lin =-| /,(p) cos n<pd<p 

- J B n B 1 »+5_„Sr n =6 1 . n =i f^^) sin n*fc 
and 

1 P 2 "" 1 f ar 

Similarly from equation (12), 

1 f 21 

A n R 2 n +A-„R 3 "=a 3ta =— ji(<p) cos n^dp 

*Kj o 

i r 2 ' 

— B n B a n +5_S s - B =6^=- fi(<p) sin 



(14) 



and 



If 21 If 2 



■02 



(15) 



•^-s xbvn ytaxn > and -D_„ THTi; /TO N n 

Let 



(16) 



i?2 

where 

T =<ri— <r,(=a / +/3') 
Substituting by means of equation (16) in equation 
(11), it is seen that 

<?l(y)=-Bo+Zj^ gnr_g-nr H e nr_ e -nr J 009 

^/di n e" r — a 3> » . ai..e _Br — 02,„\ . 
+S \ '^_ e - n / + > f _ e - Br " jsin 

?i(p)=Bo+s[-»i..(^±^) 

+M^F^)] C0S7 ^ 
+2[a,.»(^^)-a2,,( in 7j i = 5 - r )]sin 

or also 

CD 

gi(?0=i?o+S( — coth 71T+&J,, csch nr) cos 



+S( a i.» coth nr ~ °i,n csch nr) sin nip (17) 
i 

Similarly, by substitution of equation (16) in equation 
(13), 

CO 

f72(¥>)=#o+S(— &i,n csch rn-Hi* coth nr) cos ntp 
i 

CO 

+2j( a i,» csch nr— eta,, coth nr) sin ftp (18) 



The equality a 1 , 0 =aa,o=-4o is a condition of uni- 
formity that is necessary since h(z) is a regular analytic 
function in the ring region. There is, in addition, an 
arbitrary element in equations (10) to (13) (which 
may be chosen in a number of ways) viz, there is at 
our disposal the choice of the point in the z plane that 
shall correspond to, say w=c. This choice, which 
will be introduced at a later point (p. 14), is z=c 
when w—c, and will fix the constants A 0 and B 0 in 
terms of Ri and R 3 and the remaining coefficients. 

From the first parts of equations (14) and (15), 

Denoting the variable in equation (17) by <p' instead of <p and substituting by means of equations (14) and 
(15), it appears that 

CD 

ffi(p / )=-Bo+~^ f /i («»)(— sin n<p cos w'+cos n<p sin n<p') coth nrdip 
i 

+~^C /s(*0(sin n<p cos n<p'~ cos n<p sin nip') csch mdtp 
i 

or 

gi(y / )=-^o+^^[" J o fi(<p) sin n(<p— </) csch nrdip— J^JiM ™ ^Op - ¥>') coth nrdtpj (19) 
In a similar manner, 

ff2(y / )=-Bo+^;^j~ /2(p) sin »(?>— «p') coth nrdip— /i(?0 sin ft(?>— c/) csch 7iTd(/T| (20) 
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CO 

a) sin n{<p—<p') coth 7ir 
i 

6) 2j sin n(tp— <p') csch nr 



that occur in equations (19) and (20) may be evaluated 
in terms of elliptic functions. Consider the expansion 
for £ (u) (reference 10, p. 403) 



■ku . 27rV \ g 2 " . mm 

2t0i t0, / ll—Cr* 8111 Cd! 



In order not to confuse the periods occurring here with 
those of the preceding section, the real period is 
denoted by 2w 1 =27r, and the imaginary period by 
2w 2 =2iT, 

q=e -=e-r=^ 

Then 



i 



^fSin UU 



=^sm™ + ^±^-l) 



=y, coth nr sin nu 
i 



sin nu 



(21) 



Consider the expansion for f(u+« 2 ) (reference 10, 
p. 426) 6 

e'(«) 



f(«+« 2 )-^-7,'=Z(«)= 



e(«) 



27rV^ a" . Tia-ii 

= — > , * „_ sm 

tor/ il—q 2 * wi 



For 6)!=^ this expression becomes 

CO 

Z(u)=^ T 



2e _nr 

— z^sinnu 



=y~! csch tit sin nix 
i 



(22) 



Then replacing « by equations (19) and (20) 

become 



1 For reference, note the definition of the 6 function (cf. references 11 and 12) 
e(u)-l-2j cos "+2^ cos ^-25» cos ^+ . . . 

The n function Is defined In equation (2.6) of the preceding section. Some writers 
on elliptic fractions use So and 0i to denote the 6 and H fractions. 



"JO 

i C 2r 

— -J fi(<p)Zi(<p—<p')d<p 

T Jo 
1 f 2r 

— -J fi(<p)Z(<p—<p')d<p 
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(23) 



(24) 



Since Z(«)=-|^ and ^i(«)=|^, there is obtained 
also by integration- by parts: 

0i (*>') =5 0 - \ J* ft W) log e { v - <p')d<p 

+1 fVte) log B.(<p-<p')d v (23') 
ff s (^')=5 0 -^ f 2 /,'^) log H.( v - V ')d<p 

7TJ0 

+;J^/i' W lo s 9 (<p~<p')<b (24') 

where the logarithm operates only on the absolute 
value of the quantities 6(p— <p') and ~B.(<p—<p'). 

In a manner similar to the foregoing procedure it is 
possible to solve for the coefficients in equations (11) 
and (13), substitute in equations (10) and (12), and 
obtain as the reciprocal relations to (23) and (24) the 
following: 

i r 2r * 

M<p')=A 0 — - gi{<p)Z{<p-<p')d<p 
« Jo 

+-J o gAdZAv-v'Wv (25) 

i P 2 * 

/ , s( ¥ >')=A — ff«(«p)Zi((p-V)<iV 

PVi(*0Z(*- (26) 

Equations (23) and (24), which essentially express 
a pair of boundary-value relations for a concentric 
ring region, permit the obtaining of the imaginary 
parts of a complex function h(z) along the boundary 
circles of a ring region from a knowledge of the real 
parts along the boundaries. These equations are 
fundamental in, a potential-theory study of ring regions ; 
they have been developed in a different manner and 
for another purpose by Henri Yillat in 1912 (reference 
13, p. 147). It will be shown shortly that equations 
(23) and (24), when generalized and regarded as 
integral equations instead of being considered as 
definite integrals, make it possible to obtain the com- 
plete correspondence which we are seeking for doubly 
connected regions. 
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The function giving the value of h(z) at any point interior to the ring region may also be expressed in terms 
of the real parts of h(z) along the boundaries d and C 2 . From equation (16), 



1 



Then equation (6) 

h(z)=a 0 +'i^ (a„3-+a_ B z- B ) 
i 

becomes 

h(z)=a 0 -- rMriMdv+l F*M<<>)Nd<p (27) 



where 



M_ ? 2 sinh tit 

^ V 2 sinh tit 
The quantities M and N may be readily expressed in 
terms of elliptic functions. Let 

Then by equation (22), 



=iZ(v 2 ) 



M _»ljll£, = e'(t) 3 ) 
^ sinhnr i0(o 2 ) 
Similarly let 

Then 

Then finally 

A(s)=a 0 -- f 2r f l (< P )Z(v 2 )d v +- r Wp)Z{v,)d v (28) 
where 
Vi=i log 

or writing 

»i=¥>— <p'+i(<r— (Ti) 

Determination of the constants A 0 and B 0 . — It will 
be recalled that the neighborhood of the point w=c 
corresponds to the region at infinity in the t and u 
planes. In order to make the correspondence of the 
w and z planes unique the following condition is put 
down. Let z=c when w=c, hence causing the region 
about z=c to correspond also to the region at infinity 
in the t and u planes. There is, however, an essential 

fact to be noted, viz, evaluated for z—c is, in gen- 
eral, different from unity, hence generally a magnifica- 
tion and rotation of the regions near w=c and z=c 
exists in the two planes. 



The conditions to be studied are 

^=re <£ , evaluated for z=c 
dz 



(29) 
(29') 



From equation (5) 



there is obtained 



(30) 



The condition (29) then corresponds to 

[h(z)] t . c =0 (31) 

And, in view of the preceding relations, equation (29') 
corresponds to 

[ Z ^] tm =re^-l=p+ig. (31') 

where it is noted that r and £ are given in terms of 
p and 2 as follows: 



£=tan~ 



By equations (6) and (9), it is found that equation (31) 
separates into 

^o+S(^»c"+^--c-)=0 



5 0 +Z;(5,c"+ J B_ n c-»)=0 
i 

Also equation (31') becomes 
i 

S7i(B„c B -5_ B c- B )=2 



(32) 
(33) 

(34) 
(35) 



These equations may also be expressed in other 
forms. For example, from equation (28), since z=c 
corresponds to <r=0, <p'=0, we have that 

+- f M<p)Zi(vi)d<p 

where 
and 

Vi=<p — icy i 
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Employing equation (22), this separates into 

CD 

. l\~Vsinh 7i<r 2 f 2 *", , . 7 
1 

oo 

+S IeS D M cos 7i ^ =0 (32,) 

i 

CO 

^-SSaSl 2 ^ ^ **** 

1 
1 

In these equations/! and/ 2 may each be altered by the 
addition of a constant without altering the values of 
the integrals. Hence, if fi(<p)— A 0 , j 3 {<p)—Ao are 
known (i. e., only the variational parts of / t and/ 2 are 
known) and if there are also given or known the quan- 
tities Ri=ce'i, R 3 =ce?i, equations (32') and (33') de- 
termine directly the constants A 0 and B 0 so that con- 
dition (31) is satisfied. 

Since /ii and /* 2 differ from / x and / 2 by constants 
(cf. equations (7) and (8)) they may replace/i and/ 2 in 
equations (32') and (33')- Also, by equations (14) 
and (15) it is recalled that 

1 C 2 ' 1 C 2r 

l f 2 * 

orthatr=(Ti — a>= 2v J — ^ d( ^' •^■ eiice ' ^ there 
are given the functions pi and hi (therefore r is known), 
equation (32') determines the individual quantities 
tri and <r 3 Q.. e., the radii Ri and i? 2 ). Equation (33') 
then again defines the value of the constant B 0 . 

By the use of equations (15) and (16), the conditions 
(34) and (35) may also be written in other forms. 
Thus 

ca 

V t= *1j sinhnr Jo U<p) C ° S 
i 

x 

00 

1\~V n sinh «<r 2 f 2r , , . . , 

2=^2/ sinhnr J 0 ^ Sm 
i 

~S Z? ^nF^? i Jo /a M ^ n<pdip (35,) 
i 

Further study of equations (23) and (24). — It has 
already been mentioned that there are two points of 
view from which the simultaneous equations (23) and 
(24) may be studied. In one, the equations are re- 



garded as definite integral evaluations and the func- 
tions fiM and / 2 (p) are known as functions of the 
variable <p. In the other, the equations are regarded 
as integral equations and the functions are known in 
terms of 8, not <p. In the next few paragraphs the 
definite-integral viewpoint will first be employed and 
it will be shown how it may be used to develop biplane 
arrangements in an artificial or indirect manner. The 
results obtained by this method will also be of some 
interest and value when the more direct integral- 
equation point of view is investigated in the subsequent 
section. 

Families of biplane arrangements. — When/i(?>) and 
ji{<p) are known functions, the evaluation of equations 
(23) and (24) determine the "conjugate" functions 
gi(<p) and g»(<p). It may be observed that the Fourier 
series expansions of fifo) and gi(<p) and of / 2 (p) and 
g*(<p) are related by the peculiar interchange of co- 
efficients as seen in equations (10) to (13). The 
existence of the integrals in equations (23) and (24) 
requires only that fi(v) and / 2 (p) be piecewise con- 
tinuous and differentiable, and have no poles of order 
equal to or greater than one. In this paper, however, 
the only interest is in continuous, single-valued func- 
tions/! and/ 2 , of period 2ir, and satisfying the condi- 
tions of uniformity (cf. paragraph following equation 
(15)), such functions as may always be associated with 
the conformal transformation of doubly connected 
regions bounded by continuous closed contours for a 
proper choice of coordinates. 6 

When the functions gi{cp) and g 3 (<p) are known, the 
correspondence of 9 and <p is imm ediately known along 
the boundary contours since g(<p)=d—(p. Also, the 
functions /i (^) and/ 2 (¥>) together with the constants 
<ri and <r 2 determine 7 the functions ni(<p) and nz(<p). 
The quantities pi and ju 2 expressed as functions of 
0(=— X) then permit the defining of two contours in 
the u plane, the external region of which is in one-to-one 
correspondence with the ring region. Some specific 
examples will shortly be given. With an insight 
gained by experience, the functions ni(<p) and /x 2 (p) 
may be so chosen that certain desired classes of prac- 
tical biplane arrangements may be obtained. It may 
be remarked here that once there is obtained a definite 
biplane arrangement by means of this process, the 
problem may immediately be considered reversed and 
thus insight is obtained into the solution of the asso- 
ciated integral equation. This notion will be later 
examined; in this section, some illustrations of the 

afore-mentioned process are briefly presented. 

* 

• There b an additional condition on ffiO) or mW) necessary in order that the 
contours shall be free of double points (cf. reference 1, p. 10), vir, 

' It Is understood that the minor equation (31) is to be satisOed. 
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or finally, 

Mi ($0=1.4708— 0.1 sin <p=- 



By reference to the Fourier series developments for 
fi(<p) and ji{<p), equations (10) and (12), it may be 
observed that a particularly simple example is the 
following: 

/i(?0=Mi(¥>)— <ri=A 0 — 0.1 sin <p 

=A 0 — (BiBi—B-iBr 1 ) sin <p (36a) 

/a(v)=M2((p)— 0-2=^0+0.1 sin <p 

=A 0 -(B 1 R 2 -B_ 1 R 3 - 1 ) sin <p (36b) 

Also let 0-1=!— 0.1=1.4708, 811(1 let °" s=— i- 4708 - 
(Hence t=o-i— <t s =2(1.4708) and 2w a =2iT=5.8832i.) 

y 




(7-0°) 



9ii.<p)-- 



— 0.04851+0.11114 cos v =g 3 (<p) 



It is now possible to define the variable 0=tp+g(<p) 
along each contour: 

0i= — Xi=<?+0i(<p) 

0a= — Xj=*>+f72 (*?) 

In addition to the choice of /i and/ 2 , the line seg- 
ments, or chords, and their relative positions (deter- 



(r - 90°) 



lb) 



(a) 




o 














(r-30°) 



(c) (d) 

Fiquee 6.— Biplane arrangements defined by special choice of fi and fi and other parameters. (See equations (36a) and (38b).) 



By equations (11) and (13), there may be written for 
#1 (?>) and g 2 (<p) 

g 1 (<p)=(6-<p) 1 =B 0 +(B 1 R 1 +B- 1 Rr 1 ) cos v 

g 3 (<p)=(6—< P ) 2 =B 0 +(B 1 R a +B- 1 R 1 - 1 ) cos v 

The conditions on A 0 and I? 0 , equations (32) and (33), 

0 1 T 

give Ao=0; Bi=B^=-^ csch^ (cf. equation (16)), 
hence B 0 = — 2Bi. (Also, note that for thin example, 
equations (34) and (35) give ;p=2=0, i. e., j^gfj 
=re'£=l). Then, 

9i(fp)= — 0.1 cschg+O.l cothgCos <p=g2(<p) 



mined by a, /3, y, and 5) about which the biplane con- 
tours will be generated may be chosen. (See fig. 5.) 



Choose 5=0 and a=/3=^ (hence 2u'=2iri). 



The rec- 



tangular Cartesian coordinates of the biplane contours 
are now given by equation (2.29). Figure 6 shows the 
arrangements obtained in this numerical case for a 
number of values of the angle of stagger y. 

In the elementary example just described the profiles 
for zero stagger are mirror images (fig. 6 (a)). A simple 
numerical case for which this is not true is given by 

Mv)=fii(.<p)-<ri=A 0 -0.1 sin (?+30°) (37a) 

M<p)=^(<p)— 0-2=^0+0.1 sin (<p— 30°) (37b) 
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Here, also, choose a=P—~i <r l =—o- 2 =lA708. Then 

gi(,<p)=(6— ¥ >) I =5o+0.09625 cos 0.0450 sin <p 

gt(<p)=(d— f >) a =5o+0.09625 cos p+ 0.0450 sin <p 

The constants A 0 and B 0 determined, by equations (32) 
and (33) are 

^io=0.0216, 5 0 =— 0.0420 



(Again, in this example jj^-J 



, re 'f=l- ; ) 



Forming di and 0 2 , the rectangular coordinates of the 
contours are obtained as was shown in the preceding 



Indeed, if the process is considered in the light of a 
boundary-value problem of the concentric ring region, 
it is seen that it is sufficiently general to yield any 
biplane arrangement (more generally, any doubly 
connected region). 

Equations (23) and (24) as integral equations. — It 
is desirable at this point to introduce the following 
notation. Frequent use will be made of subscripts. 
The first subscript will usually be 1 or 2 and will indi- 
cate that the designated quantity is to be evaluated at 
the boundary C\ or-Ci, respectively. (ox-also-jBi-and^,, 
respectively). A second subscript will sometimes be 
employed to denote the variable in terms of which the 
quantity is expressed. Thus ah.s represents the quan- 




(r--l5°) 



(r - -30') 



o 


o 












(7-15°) 




(7- 30") 




W ( e ) m 

Fkiobe 7.— Biplane arrangements defined by »pecial choice of ft and/i and other parameters. (See equations (37a) and (37b).) 



example. Figure 7 shows some arrangements for 
various values of the angle of stagger 7. Figure 7 (f ) 
shows the arrangement obtained (for 7=0, cf. fig. 
7(a)) when an angle of decalage 5/2=— 3° is further 
introduced in this numerical example. 

In this manner, by employing appropriate values for 
fx(<p) and/ 2 (pO (or g t (<p) and g s (<p)) and for the other 
parameters involved, it is possible readily to develop 
arrangements of a great variety of contour shapes, 
gap/chord values, chord ratios, stagger, and decalage. 

130092—37 8 



tity /x evaluated around C x (or Bi) expressed as a func- 
tion of 8; also, /i 2it , denotes the quantity n evaluated 
on C% expressed as a function of <p. 
It is recalled that by equations (7) and (8) 

/l.p—Mi,* - ""l 
fi.r—Pi.r — ff * 

orfi and hi, and/ 2 and 11% differ by constants. Then 
in the integrands of equations (23) and (24) the func- 
tions/! and/ 2 may be replaced by jtii and /i 2 , since the 
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additive constants do not contribute to the integrals. 
If a definite, biplane arrangement is preassigned or, 
what is essentially the same, if there is given a definite 
nearly circular annular region (in a id plane; w=ce" +ia ), 
it is the functions nu and ^a.o that may be considered 
directly defined or known. 8 Thus, from an initial 
knowledge of mi . e and jt 2 , e and with the aid of equations 
(23) and (24), it is desired to obtain a knowledge of 
the functions pi, p and n 3 , f . From this point of view 
the expressions (23) and (24) then represent a pair of 
simultaneous integral equations, whose process of 
solution is more intricate than that involved in the 
process of evaluating the definite integrals. The 
problem may be restated more precisely: 

Given two functions jii.9 and /i 2 ,« that define two 
continuous contours of an annular region with respect 
to an origin contained by both contours. Then two 
pairs of functions iii t<f and g it9 , and /n 2 ,„ and gi, v are 
to be obtained such that they are interrelated in the 
manner shown by the interchange of coefficients in 
the Fourier expansions, equations (10) to (13), which 
also satisfy certain local specified conditions at z=c 
(equation (31)), and for which the relations g t = (0— ¥>)i, 
ga=(d—<p)2, which permit an interchange of the argu- 
ments 6 and <p on each boundary, are consistent with 
the given functions jti,« and that is to say, when 
Pi itf is expressed as a function of 0 by means of the 
function g ll9 there results the original function p u t; 
in symbols 

Hi, ? = 9 

when 8=<p+g 1 , f or <p=0—g lw e 
and similarly 

when 6=<p+g 2wV or <p=d— g 2 ,e 

The process employed in this paper to obtain the 
desired solution of the simultaneous integral equations 
is one of successive approximations or iteration. 
The degree of convergence of most methods of suc- 
cessive approximations usually depends on how good 
the initial approximation is. In this regard it is 
rather fortunate that the contours of practically any 
biplane arrangement transform, under the transfor- 
mations already developed (with proper choice of 
coordinates), into nearly circular contours in the to 
plane. The nearness of these boundary contours to 
circular contours is very significant and enables the 
initial approximations to be so chosen that the process 
converges ordinarily with great rapidity, one step in 
the process being sufficient for most practical purposes. 

Outline of the method of successive approxima- 
tions. — The various steps in the process of successive 
approximations will be written down schematically. 
(Allowing for a difference of notation the process is 
essentially similar to that employed in reference 1.) 



• The process of forming /ij , and assuming as known oniy the rectangular 
coordinates of the biplane contours, is outlined In section 4, p. 68. 



An extension in the use of subsequent notation must 
first be noted. The symbol represents, as men- 
tioned previously, the function / evaluated on & and 
expressed as a function of 0. The symbol /i.t,« shall 
now be employed to denote the value of /i as given 
in the ktk step in the process of successive approxi- 
mations, expressed as a function of 0. Thus the 
symbol g2.t. n _ 3 denotes the value of g 3 given by the 
fourth step in the process, expressed as a function of 
<p 2 as given by the third step in the process. 

We start with the two functions jti.j and n 3 .e that 
define the contours B x and S 2 completely. Employ- 
ing 0 instead of p in the integrals, 9 equations (32') 
define the constants f U and a- SA and equation (33') 
determines the constant B 0A . The simultaneous 
equations (23) and (24) then determine completely 
the functions 51.1.9 and <7 2 .i,9 as follows: 

9i.i=e-gi.i.i 1 (38a) 

¥>2.1 = 0— 02.1.9 J 

The values hi, 91A and m, n i may now be defined 
and these functions may be considered as known. 
Employing <pi A and tp%.i as variables instead of <p, equa- 
tions (32') and (33') determine tr 2 . 2 , and B 0 . 2 . 
Equations (23) and (24) determine then the functions 
g lm2> n i and g».3, r i V which may be expressed as func- 
tions of 0 in view of equation (38a). The variables 
^i.a and <p3.t are now given by 

<Pi.*=e-gi.*. 9] (38b) 

VS. 2=0— 02. 2, 9J 

The functions jm, ^ 2 and p % n . s are now determined 
and the process may be continued as outlined for 
H, t n j and /i 8 , n r It is noteworthy that this process 
converges, in practice, with extreme rapidity, that is 
to say, the functions ni, n ^ and na, n . k approach iden- 
tity with Hi.n.jt+i ^ to.n.tH for 3111011 values of k. 
Experience has shown (cf. also reference 1) that in 
ordinary cases one, or at most two, steps in the process 
are sufficient for great accuracy. It must be noted, 
however, that a completely rigorous discussion of the 
convergence process is lacking. 

In order to illustrate the method, consider the 
biplane arrangements of figure 7 defined by the func- 
tions in equations (37a) and (37b). Forgetting for 
the moment that the various functions are known, 
assume only pi,s and /t2,9 to be given and attempt to 
obtain p u r and /is. 9 - Figure 8 shows the various 
functions described. It is seen that the set of curves 
ffi.m.i. 92.n,i. and /t 2 . W-1 (obtained by a 20- 

point numerical process similar to that sketched in 

> Denoting the initial approximation by a rero subscript, observe that the Initial 
appradmationemployedhereIs?>i.t— 9— ffi.iand (oi.o^s— where ffu)—0j.o=O. Moro 
generally, the initial transformation may be better defined where <J\,% and flj.o are 
arbitrary functions, so chosen that they are better approximations to the Dual 
solutions ffi and gj. Then employing <n.t and as variables Instead of 9 the func- 
tions wxfu and pi.i. n j may be defined from which nij and n.i are determined 
and the process continued as outlined. (Of. reference 1, p. 13.) 
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the appendix of reference 1) axe completely coincident 
with the known solutions g iir , g 2ir! ni fV , and n 3tf . A 
further application of the process can cause no fur- 
ther noticeable change. 




0 30 60 30 120 150 180 210 240 270 300 330 360 
8,p , degrees 

FiauBE 8.— Illustrating the process of successive approximations. 

For future reference, the derivative expression ^ 

evaluated at the boundary circles will be required. 
From equation (30) 



dw w, 
~3z 



On Oi 



Then 



_wf _ dh{z) \ 



Observing that d<p=d(8—g) this may also be expressed 



dvT\ 



(39') 



For the boundary C 2 , the subscript 1 is replaced by 2. 

4. GENERAL MONOPLANE WING SECTION THEORY — 
A DEGENERATE CASE OF THE BIPLANE ANALYSIS 

It may be of some interest to discuss in this section 
briefly the case of the single airfoil section considered 
as a special case of the biplane analysis. A biplane 
arrangement in a two-dimensional field of flow corre- 
sponds mathematically to a doubly connected region. 
The degenerate case in which one of the biplane con- 



tours reduces to a (regular) point leads to the mono- 
plane-wing profile the external region of which is 
simply connected. 10 It should therefore be possible 
to obtain the complete theory of monoplane wing 
sections in potential flow as limiting values of the 
formulas already developed for biplane wing sections, 
as will be outlined in the present section. Conversely, 
the monoplane case is of further significance in that it 
permits a more complete understanding of the biplane 
analysis and, too, considerably simplifies the practical 
evaluations. The numerical process employed in 
the monoplane case, it will appear, needs to be modified 
in a relatively minor way to yield the results for or- 
dinary biplane combinations. Indeed, according to 
a method outlined by Koebe (reference 2), it is known 
that even the more general problem of transforming a 
multiply connected region (multiplane problem) into 
a region bounded by circles may be resolved by a 
process of successive approximations employing only 
the separate cases for the simply connected regions. 
The details of this problem remain for a future investi- 
gation; it is mentioned here merely as a further possible 
application of the degenerate case. 

It may first be observed that as one of the circles 
in the t plane reduces to a point the value of a or f) 
(according as the upper or lower circle degenerates) 
becomes infinite. Hence in the degenerate case, the 
elliptic functions introduced in the analysis become 
circular functions having a real period 2w=2w and an 
imaginary period 2w', which is infinite. From equa- 
tion (2.23) it is clear that &/=•£<» corresponds to g=0 
and that 



Lim A(s)=^ cot|+| cot^ 



g=0 



(1) 



Let jS=co, then employing this equation and relation 
(1.8) it is found that equation (2.12) becomes 



Limu(s) e" /a =Z7 I =T 1 +Sr 
2=0 Jl 



(2) 



where the origin in the Ui plane is referred to the mid- 
point of the line segment and where 

r 1 =e«w«m(t_i c coth a) 

a=c csch a 

It is seen that the T x plane corresponds simply to the 
t plane translated to a new origin (Ox in fig. 1) and 

rotated by an angle 7' =7+^* Equation (2) is noth- 
ing more than the well-known transformation leading 
to Joukowsky airfoils. The line segment in t.hia case 
(|3=co) is equal to twice the diameter of the circle in 
the T x plane, i. e., Ci=4a. 

Strictly speaking, the point at Infinity represents another boundary and the/tow 
regions are respectively, triply connected and doubly connected. This tact corre- 
sponds to the circumstance that in the flow formula for the biplane case there may 
be two arbitrary circulations specified, and In the monoplane flow formula, one. 
The flow formula for a simply connected region without singular points cannot be 
multiply valued, 1. e., cannot possess an arbitrary circulation. 
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If a instead of (3, is allowed to approach, infinity 
there is obtained 

6 2 

Lim u(js) = U a = T 3 +7p- 

where the origin in the Z7 2 plane is referred to the mid- 
point of the line segment, and where 

T 3 =e»(t+ic coth 0) 

b=c csch p 

The chord length c a =4&. 

In what immediately follows we show how to express 
X, /i in the degenerate case in terms of the coordinates 
x, y of the airfoil section. These results will be of 
value in the later determination of X, ix from x, y for 
the biplane case. 

It has been seen that in the degenerate case (j3=q>) 
the origin of coordinates may be referred to the mid- 
point of the chord and 

where the rectangular coordinates in the TJ X plane are 
(x u yi) i. e., 

U^xi+iyt 
Tt=err'(t—ic coth a) 
a=c csch a and y'=7+ 5/2 
Let Ti be written in the form 

!T 1 =oe*i + "i (3) 

Then 

U l =x 1 +iy l =2a cosh (fi+idi) (4) 

where 

x 1 —2a cosh cos 0i 
y!=2a sinh ^ sin 0i 
And upon inversion (cf. reference 1) 



2 sin 2 0,=p+yp*+| 



where 



2sinh s ^=- ? +y?+| 



(5) 



By equation (1.8) 

f=— c cot | where s=X+in 

Hence, 

X+i M =-2 cot-'^=-2 cot-KTie-n'+ic coth a)^ 
or also 

\+tn=-2 cot" 1 (Z+im) (6) 
where 

l=-e*i cos (0i—7') 
c 



m=^(e*i sin (fii— T^+cosh a) 



It is known that 



cot HH-tm^-^log^pj 



-im— t 



(P+m*-l) 2 +4P , 1 x _! 2Z_ 
- * 4 wg [fJ+(m _ 1)2]S i- 2 tan p+m! _ 1 

So that finally equation (6) separates into 



X=— tan" 



, 21 
Z 2 +m s -l 



1, (P+m s -l) s +4P| 
"=2 log [P+(m-l)l* 



(7) 



These relations express X and n in the monoplane case 
in terms of B t and and hence by (5) also in terms 
of Xi and y t . 

Similarly, in order to obtain X and jtt for the degener- 
ate case in which <*=<», replace a by 6 (=c csch |3) 
and let 

T 3 =e" (t+ic coth 0) 

=&e*+«a (3') 

Then finally X and n are given by equation (7) in whioh 



I— Z 2 =-e*» cos (0 a — 7) 
c 



m~m i =^(t* 1 sin (0 a — 7)— cosh 



(60 



Inversion of equation (2.12). — It is quite evident 
that the direct inversion of the elliptic transcendental 
equation (2.12) (p. 52), if at nil possible, would be 
very laborious. However, an indirect method which 
employs the results of the degenerate cases and which 
performs this inversion readily, to any degree of 
approximation, will now bo outlined. 

Thus, let a definite biplane arrangement be given 
(fig. 5(a)). The chords or line segments c x and c a may 
first be chosen in the following convenient manner. 
Let the chord c x be defined by the line F\F\, where Fi 
is the midpoint of the distance between the leading 
edge and the center of curvature of the leading edge 
and Fi is the midpoint of the distance between the 
trailing edge and the center of curvature of the trailing 
edge. It is observed at this point that the only theoret- 
ical restriction upon the choice of the chords is that the 
singular points (the end points of the line segments) 
be within the contours, or, at most, on the boundaries 
themselves. The above-mentioned choice is one 
merely of convenience, the object in view being the 
defining of a smooth (X, p) relationship. (In reference 
1, a similar situation is described in detail.) The 
above-outlined procedure may be also applied to 
determine the chord c% of the lower contour of the 
biplane cellule. 11 

it In the event that the chords so determined axe almost, bat not quite, parallel 
it Is ot some advantage numerically to vary from the foregoing procedure suf- 
ficiently to cause the chords to become exactly parallel and to maintain approxi- 
mately the foci F and J". Small variations from the choice of chords outlined are of 
minor importance and will not affect the smoothness of the (X, ft) relationship. 
The desirability of maintaining the chords exactly parallel, If they are reasonably 
parallel to start with. Is due to the circumstance that elllptlo funotlons of the sec- 
ond kind are then avoided. 
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The chords Ci and c 3 having been conveniently chosen, 
it is possible to determine uniquely by means of the 
charts outlined previously (fig. 4), or indirectly from 
the equations themselves, the values of the constants 
a, /3, 7, and 5. The quantity c may be regarded 
throughout as a convenient unit reference length. 

Equation (7), it will be recalled, for these values of 
a, /9, 7, and 5, determines the values of X and y (in 
terms of the rectangular coordinates of the profile 
sections), which in the degenerate cases correspond to 
profiles geometrically similar to the upper or lower 
profiles of the given biplane cellule, the only differences 
being that the chords are la=ic csch a and 4b= 
4c csch /3, respectively. When, however, tho values of X 
and y thus determined are inserted in equation (2.29) 
(employing tho proper periods 2co=2r, 2co'~ 
2i (a+/3) ) there is obtained a biplane arrangomont 
which has, necessarily, the required chords and position 
(i. e., tho proper skeleton) and around this skeleton 
has generated contour shapes A x and B x which, in 
ordinary cases, are almost identical with the given 
original contours A Q and B 0 . Thus, the use of the 
values (X, y) of the degenerate cases in the biplane 
analysis is equivalent to a replacement of the original 
biplane arrangement by a new arrangement defined by 
tho contours A x and B x . The differences between A 0 
and Ai f and B 0 and Bi are remarkably small in prac- 
tice. Mathematically the foregoing procedure repre- 
sents, however, only a first, although important, step 
in a process of successive approximations which we 
outline as follows. 

Consider only A 0 and A x . The contour A x defines, 
by means of equation (7), a new degenerate (X, y) 
relation. The differences between this new (X, y) 
relation of A x and the original degenerate (X, y) relation 
of A 0 is a proper criterion of the differences between 
tho contours A 0 and A x themselves since, if these (X, y) 
functions coincide, the contours must coincide. Hence, 
by a shifting process similar to that commonly em- 
ployed in methods of successive approximations, the 
first approximation to the desired (X, y) relation of A 0 
(this first approximation has boon here considered 
to be the degenerate (X, y) relation itself of A 0 ) may 
be corrected by those differences to give a second and 
better approximation. The process may be repeated k 
times, if necessary, until the degenerate (X, y) relation 
of A t coincides with tho degenerate (X, y) relation of 
A 0 , hence A t coincides with A 0 , and therefore the 
actual (X, y) relation that defines the contour A t 
itself in the biplane case (i. e., by equation (2.29)) is 
the desired (X, y) relation of A 0 . 

Singular points in the monoplane case. — The mono- 
plane case may also be useful in obtaining the singular 
pointR of the biplane transformation to a good firsts 
approximation. From equation (2): 



U 1 =T l + 



the singular points are given by 
or by 

Since by equation (3) 

it is evident that the singular points correspond to 

(fi, B 1 )=(0, 0) and (0, tt), respectively. On replacing 

4>i and 0, in equations (6) and (7) by these values it 

is seen that the singular points (X=X,) of the upper 

chord (y=a) are given by 

. . sinh a cos y' 

—tan X,= q— ; -. ^ — 7 (8) 

1± cosh asm y' K J 

Similarly for the chord of the lower profile (ji=— p) 

we get 

. . sinh jS cos y ... 

tea X,=i-; , a i 1 — (8') 

1 ± cosh 0 sm 7 K ° ' 

It may also be useful to note that to a first approxi- 
mation (the approximation being better the greater 
a+j3) tho chords of a biplane cellule are given by 
Ci=4c csch a, c 3 =4lC csch j8. Hence the chord ratio is 
approximately 

Ci_ sinh 
c s sinh a 

It has been shown thus far that the degenerate mono- 
plane case may be of value in the determination of 
y and X for a biplane cellule. The limiting forms of 
the simultaneous integral equations will now be briefly 
discussed. 

Forms of the integral equations (3.23) and (3.24) 
for the monoplane case. — Consider the simultaneous 
equations (23) and (24) of section 3, wliich define 
the distortion of two contours from two circles. 
Let (3'=— <rj=co, i. e ., B 2 =0, the interior circle of 
figure 5 (e) degenerates to a point. It is seen then 
that the Laurent series, equation (3.6), becomes a 
one-way ascending power series and o_ B =A_„+iB_ B =0. 
We have also since T=a'+j3'=co that equations (3.21) 
and (3.22) give 

Z( p -V)=0, Z,{v-S)=\ cot 

Equations (3.23) and (3.24) then reduce to a single 
equation 

gi{<P , )=B 0 -±J*f l (<p) cot ^^<h (9) 

or also 



gi(v')=B 0 +l f*M v ) log 



sm 



d<p (9') 



Equation (9), considered as an integral equation, 
enables the transforming of tho contour of a simply 
connected region into a circle. The process of iteration 
outlined in the preceding section is directly applicable 
in this simpler case. (Cf. reference 1.) The constants 
A 0 , B 0 , and or,=a' may be determined, as before, from 
equations (3.32) and (3.33). 
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In the event that a'=<°, the outer circle of figure 
5(e) becomes infinitely large, i. e., J?i=t». Then the 
Laurent series (3.6) becomes a descending power 
series and a n =A n +iB n =0. With /i=<h=0 and 
T =oo -we find that equations (3.23) and (3.24) reduce to 



or also 



92(.<p')=B 0 +^ c \ f 2 (<p) cot— 2~ d<? 
so 

1 P 2 * 

g*(<p')=B a --J a / 2 («p) log 



(10) 



d<p (10') 



The functions gi(tp) and g^tp), determined by equa- 
tions (9) and (10) by separate treatment of the two 
monoplane cases, may when known be employed as a 
convenient initial approximation (cf. footnote 9) in 
the more general case. 

The well-known flow function for a single circular 
cylinder may be brought into combination with the 
results of_this section to yield the flow about an arbi- 



spond to the contours of the biplane arrangement and 
are not to be confused with the contours Ki and K 2 of 
the t plane (fig. 5 (d)) which, in general, are not circles. 
The primes will be retained to denote this difference. 
The relation between the z and f planes is (cf. equa- 
tion (1.4)) 



t'=xc 



z+c 
'z—c 



Also by the relation 



t'=—C COt 7T 



(1) 



(2) 



the region external to the circles Ki and K 3 ' is mapped 
into a rectangular region in the s' plane bounded by 
the lines £/ and W (cf. equation (1.8) and fig. 9). 
For later reference the relation between s' and z is also 
noted here: 

z=ce- u ' (3) 
and since z~c^ +i ' p and 8'=\'-\-ip' it is clear that 



] n z ' p/ape 




r 



r. 



s' plane 



IT*' 



n 




\ V plane 


1 ft 


<?. 


"~~\-c csch d' 


t 




Vm 


o 











(a) 



r, 



(b> (c) 
Fiouke 9.— HIostraHns the flow regions In the t, f, and f planes. 




c csch /3' 



trary single airfoil as has been already obtained in 
reference 3. We proceed at once to the more general 
case to introduce the flow function for two circular 
cylinders into the biplane analysis. 

5. POTENTIAL FLOW ABOUT THE BIPLANE CONTOURS 

Potential flow around two circles. — It has thus far 
been shown how the two contours of a biplane arrange- 
ment, in a u plane, may be transformed into two con- 
centric circles in the z plane. It is desirable in what 
follows to transform the concentric circles Cy and Cj of 
the z plane into the coaxial ones ET/ and K 3 ' in a t' plane 
(fig. 9). The circles Ki and K 3 ' will thus also corre- 



'=-<p\ 
'=* ) 



(4) 



Hence the boundary lines li and W are given respec- 
tively by 

n'=a'=<r 1 , n'=— 0' = <r 2 

where 

«'>0, £'>0 

(It will also be recalled from section 3 (cf . fine following 
equation (3.16)) that a\— <r 3 =a'+f}'=T.) 

In a noteworthy paper (reference 4) Lagally has 
given the complex flow potential for uniform flow past 
two circles. His formula makes use of the interme- 
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diate s' plane and the complex flow potential in the s' 
plane, is given by 12 

+2tc«c O [f(s0-f(s , +2i/30] (5) 

In this equation r= — (r^Tj) where Ti and IV are 
the individual arbitrary circulations about Ki and K% , 
respectively (see fig. 9 ; the positive direction of T 1 or r S) 
which is opposite to that of T, is chosen such that as 
one traverses the contour in this direction the flow 
region is on the right). The symbol r' denotes the 
"countercirculation" T'=Ti—T 2 . The velocity at 
infinity is u m +iv m . The periods of the elliptic func- 
tions are, as in section 3, 2wi=27r, 2oj s =2i(ci' +f}')=2ir. 

Equation (5) shows directly that the whole flow is 
built up by linear superposition of four separate flows, 
for each of which one of the quantities r, I", u m , and 
v m is different from zero. 

In the first partial flow with the factor r', both cir- 
culations are equal and opposite about the two circles, 

r' 

Ti——T 2 ~-^- The streamlines are the circles of the 

coaxial pencil. The velocity at infinity is zero. In 
the second partial flow with the factor r, the two cir- 

r 

culations are equal. ri=r 3 =— The velocity at 

infinity is again zero. The third partial flow with the 
factor u m is a translatory flow without circulation, 
normal to the line of centers of the two circles. The 
fourth partial flow with the factor v m is a translatory 
flow without circulation in the direction of the line of 
centers of the two circles. 

The flow at the surface of the two circles, and the 
determination of r and I". — It is convenient to define 
the velocity at infinity, whose rectangular components 
are u m and v m) in another manner, by introducing a 
magnitude V e and an angle of attack a e (fig. 9 (c)). 
Let 

Ua== — Vc cos a c \ 

v m = V e sin «J W 

" Lagnlly's formula as presented In reference 4 employs elliptic functions of periods 
at and ~ut which are conjugate to the periods <n and <n employed here, I. 8, m=u^/f, 
M-fui. The formula of reference 4 Is brought Into equation (S) by putting Lagally's 
symbols iZ, u a , v a equal to s", »„, -«„, respectively. 



or, in a single expression, 

u a +w><x>= — V c e _ta - 

By derivation of equation (5) with regard to s', we 
obtain the complex velocity function in the s' plane as 

W)=£-^{2iB-(s')-f(s'+2ii3')]+l-fe 
— 2cV e cos a c \j?(8 )-tp(8'+2ip')+2 Vl /Tr] 
-2icV e sin a c [p(s')-2'(s'+2i/3')] (7) 

This expression gives the velocity components TJi— iVi 
at any point of the rectangular region in the s' plane. 
It is a real quantity on each of the boundaries, p!=at ', 
m'=— jS' of the rectangle since these boundaries are 
steeamlines and the normal flow Vi vanishes. 
Let us evaluate W for each of the two cases: 

(1) 8'=\'+ia' corresponding to the boundary Z/ 

(2) s'=\'—ip f corresponding to the boundary Zj'. In 
case (1) we obtain 

-2cV c cos * c \2(\'+ia')+p(\'-i a ')+2Ti 1 /Tr\ 
-2ieV e sin cc e [p(\'+ia')-p(\'-ia')] (8) 

or 

Wl= ^r~2^ Bl ( y > a ')- 2cV < cos «A(*'f «') 

-2cV e sin acR s (y, a') (9) 

where B u R 2 , and B 3 are real quantities introduced for 
brevity in notation and defined by the following equa- 
tions (the primes are dropped for convenience): 13 

R^X, a )=i[f(X+ia)-f(X-i a )]+^-^ (10) 

1 sinli cc 00 cP m 

= -2+ cosh g -cos x - 4 ?!^. 5 ^ m <* 003 ™ X 



o The developments given here for J6, Bt, and Hi are rapidly convergent when 
T-af+P' Is large, say greater than Other expansions are possible (of. reference 
10, p. 422) and, In certain cases, more desirable. For example, when r Is small, say 
less than ^ > It Is possible by a simple transformation to Interchange the real and 
Imaginary periods of the elliptlo functions (cf. footnote 4) and obtain more rapidly 
convergent developinents. 
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1— cos X cosh a 
' (cosh a— cos X) 1 



VV mg 2 
2-1 1-2 



2 s ' 



fl s (X, «)=i[p(X+za)-j>(X-ia)] 

00 

sin X sinh a , \~\ mq lm 



(cosh, a— cos X) 1 



(ID 

; cosh ma cos mX 
(12) 

sinh ma sin mX 



where 



2 =e~ 



In case (2) there is obtained similarly 

+2cV e sin aJBtQi', F) (13) 

In order to specify the circulation r and the counter- 
circulation T', we make use of the Kutta-Joukowsky 



condition for finite velocities at the sharp trailing edges. 
Equations (9) and (13) must vanish for the particular 
values of X' that correspond to the trailing edges of the 
upper and lower contours of the biplane combination. 
Let X/ and X 2 ' be the values of X corresponding to the 
trailing edges of the two contours. 14 Then we have 



2ir 



R u —2cV e cos OeBj! 



—2cV c sin a c R sl =0 

W(W, n=^+^S»-2eV e cos ae R„ 
+2cV e sin a c R 3i =0 



(14) 



where the R's with double subscripts are constants 
defined as follows. (See equations (10)— (12).) 

2?n=B 1 (X 1 ' > a'), R tl =R 3 (W, a'), Ru=R 3 (\i', a') 

Bi»=Bi(v, n, Bm=Bi(x s / , n, B M =i^3(x 2 ^ n 



r r' 

Solving for ^ and |^in (14) 



r ^^.-f cos adRn— E M )+sin a. (R 3 i+Rm) ~\ 
2^ = - 2cV i SITfSn J 

I" T ,r cos a t (ig„i?M+.g»iBi2)-sm a e (Jg 11 g 3 »-.g 31 E 18 ) '1 
4x =2cy i BHT^ J 



In order to obtain the angle of zero lift j5„ in the 
plane of the circles, we equate 15 r=0 and solve for the 
particular value of the angle of attack a e , which is 
denoted as —fi e e., for a e =— /S c , r=0). Then we 
have at once from (15) 



tan p c = 



Rsi — R21 



(16) 



with this definition of j3 e , the total circulation may be 
expressed as 

r=— ivcKVc sin (a e +/3 £ ) (17) 

where the constant K is 

tv-__ 1 •Ks l +i2jj 
cos p e R n +R lt 

(In the limiting cases P'= a'= , cK is equal to 
c csch a' or c csch respectively, which are the radii 
of Ki or Kt, respectively.) 

Similarly, the angle y e may be defined as the angle of 
attack for which countercirculation I" vanishes (i. e., 
for a c =— 7 0 T'=0.) Then from (15) 



tan 7 C = - 



(18) 



and T' may be expressed as 



11 The values of V and Xi* may be determined as follows: Let Xi and Xj denote 
the values of X(— 0) that oorrespond to the end polnte of the chordi, L e., Xi and 
Xt, or —«i and —ft, are the smeular points of equation (2.12). Then the values of 
X'(- — corresponding to Xi and Xi are given by Xi' and Xj', or —pi and — where 
(cf. sec 3) pi-ft— ti. », and pi=>ft— jj. i r 

"The total lift Is, as in themonordajiecase,i;=pVr(cf.p.25)pernnltspenl6nrth. 



(15) 



(19) 



T' = -StcJV c sin (a,+Y«) 
where the constant J is 

j 1 -Rn-Rgj — i?3]i?u 
_ COS 7„ Rn+Ru 

Velocity at the boundary contours of the biplane 
combination. — Let the complex potential function in 
the plane containing the biplane contours (u plane) 
be £2, then the complex velocity function is 

m =v x -iv, (20) 

or 

dn 



3u 



where v x and v v are the velocity components in the 
direction of the coordinate axes in the u plane. Intro- 
ducing the intermediate planes, we have 



dQ_<Kl ds' dz dwds 
du~dlP'dlz'd^'dT'd^u, 



(21) 



It is first of importance to consider the changes 
that a velocity at infinity in the u plane undergoes 
when transformed to the f (or s') planes. (It will 
be recalled that corresponds to tf=°°, 8=0, 

w=c, z=c, s'=0, f=co.) Let —Ve~ ia i denote the 
velocity at infinity in the u plane and, as before, 
— V e e -te|j denotes the velocity at infinity in the t' 
plane. (See figs. 9 and 10; the velocity magnitudes 
are V and V e and the angles of attack are ao and a,, 
respectively.) 
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Then noting equations (2.10) and (3.29) it is found 
from (21) that 

V c =rV 



(22) 



where r and £ are determined by equation (3.31') 
(r is generally near unity, £ near zero), and y is the 
angle of stagger of the biplane chords. 

The angle of zero lift for the biplane combination 
is given by 

i. e., for «o=— /3 0 the lift vanishes. 

In order to determine the velocities at each boundary 
contour of the biplane combination, it is sufficient to 
obtain the magnitudes of the individual terms in 
equation (21) at each boundary. For the upper 
profile we have by equation (9) 

d«'\ Wl 

From equation (3) 

d*' 
Tz 



From equation (3.40) or (3.40') 
\dv 



or also 



r(i+%-)T 



'~dT 



From equation (1.7) 
And from equation (2.32) 

2^H (Pi cos y+Pi ' ^ 7)2 

+ cos7+<2i , sin 7) S J J 

(The subscript 1 in P u P/, Q u and Q, 7 denotes that 
the values of (X, p) which define the upper contour are 
to be used.) Hence, finally for the velocity p, nt any 
point (x, y) or (X, p) of the surface of the upper profile 
there is obtained 



Similarly for the velocity at each point of the surface of the lower profile 



(23) 



(24) 



It may be worth while to point out in resume 1 of 
equations (23) and (24) that these relations do contain 
the necessary parameters for the study of the potential 
flow in any doubly connected region. It is observed 
that the uniform stream velocity V and the angle of 
attack ao for the biplane cellule occur only in Wi or W 3 
(equations (9) and (13)) and are related to the velocity 
V c and angle of attack a c for the two circles in the V 
plane by equation (22). The circulation Y and the 
countercirculation I", also occur in Wi and W 3 . These 
circulations are, in general, arbitrary but may be con- 
sidered specified or fixed in tho case of biplane con- 
tours by equations (17) and (19). The parameters 
a, )3, 7, and 5 determine the position and attitude of 
the profile chords, i. e., the gap/chord, chord ratio, 
stagger, and decalage. The variables X, y. define the 
shapes of the contours, having the chosen line seg- 
ments as chords. The functions j and g furnish the 
means of transition from the two arbitrary contours 
to two circular boundaries (i. e., from 6, p. to <p, <r). 
The parameters a' and ft' determine the radii of the 
circles and are fixed by the local conditions at z=c 
corresponding to the region at infinity (equation 
(3.29)). 



The local supers tream pressure at points of the 
boundary surfaces, in terms of the dynamic pressure 
of the uniform stream, is given by 



where 



The general formulas can be somewhat simplified in 
certain special cases. For example, in the case of 
biplane contours described about parallel segments as 
chords the parameter 5=0 and, in addition, if the 
chords are equal, a=/3. 

Application may be made of the Blasius formulas for 
obtaining the total lift of the biplane combination. It 
is readily found that the total lift is, as in the mono- 
plane case given by L= pVT per unit span length, whero 
r is the total circulation, and the lift vector L is 
perpendicular to the direction of the uniform stream V. 
The total integrated moment as well as the forces on 
the separate contours may also be developed without 
regard to the local pressure formula (25) though tho 
expressions are not of particularly simple forms. In 
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the special case of biplane combinations composed of 
the framework of the line segments themselves, the 
formulas for the forces and moments will reduce some- 
what. These and further special applications are 
reserved for the future. 

As an example of the application of the formulas 
presented in this paper the pressure distribution for the 
biplane cellule shown in figure 10 (N. A. C. A. 4412 air- 




Fiouhe 10.— Biplane arrangement (N. A. O. A. Mil seotlon). 

foil section) is developed. The curves representing the 
9u 9i) Mi> and /z 2 functions' are shown in figure 11. The 
pressure distribution is given in figure 12 for values of 
the biplane combination lift coefficient: Cl=0, 0.5, 
1.0, and 1.5. The pressure distribution for the mono- 
plane case (cf. reference 14) is also presented for 
comparison. The numerical procedure is outlined 
under table I. 



Langley Memoeial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 8, 19S5. 
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(a) Pressure distribution for the upper wing section of the biplane combination In figure 10. 



-3.0 r 




(c) Pressure distribution for the monoplane airfoil section (N. A. O. A. 4112). 

Fiqubb 12. 
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EXPLANATION OF TABLE I 

Table I presents in outline the numerical procedure 
in obtaining the velocities and pressures at the 
boundaries of the biplane combination shown in figure 
10. The gap/chord=0.825, stagger/chord=0.5S0 ) 
chord/chord=l, and the decalage=0. The param- 
eters a, 0, 7, and 5 are: a=j3= 1.4436, 7=— 30°, and 
5=0. The rectangular coordinates x, y of the profiles 
are given in columns 1 and 2. (c=unity.) The 
values of X, m which correspond (i. e., satisfy equation 
(2.29)) are given in columns 3 and 4 and have been 
calculated by the method outlined in section 4. 
The constants A>=0.0467, #0=— 0.0106, r=1.025, 
{=— 0°3', a'=1.333, £'=1.303, t=2.636 and the 
angular distortion functions gi and g 2 (column 5) are 
deterrnined by the method outlined in section 3. 
Column 6 presents the angle <p=d—g where 0=— X. 
The functions jB,, R 3 , and B 3 given in col umns 7, 8, 
and 9 are determined by equations (5.10)— (5.12). 

Column 10 gives the quantity ^=[(! +3^) "^(i^)*]" 1 
for each profile obtained graphically. (See fig. 11.) 
Column 11 gives the quantity k=[(P cos 7+ 
P' sin (<2 cos y+Q' sin 7)*]"* for each profile. In 
determining the next column, we require the singular 
points of the chords as determined by equation (2.22), 
X,=91°49', X 2 =40°5S'. Hence Vl =— 93°22', <p 2 = 
— 35°14' (cf. footnote 14). The constants occurring 
in equation (5.14) are then JJ n =0.3444, R :1 =0.2Q17, 
i?s,=0.3659, £ 12 =0.9431, B 2i = — 0.4890, £, 2 =0.7117. 
Equation (5.16) gives then /3 C =34°51' and equation 



(5.17) determines X=1.020. The circulation is now 
r=— 4w*V e (1.020) sin (a e 4-34°51 / ). Equation (5.18) 
determines 7 e =38°5' and (5.19) gives J=— 0.0987. 
The countercirculation Y r =8ircV e (0.0987) sin (a e + 
38°5')- The lift coefficient may be expressed as 0 L = 

— -k where q=- F; pV 2 . The chord for each profile in 
g(cnora) I 

the example considered is equal to (2.106)c, hence 
C L = — 2tt(0.985) sin (a„+34°51'). (In the mono- 
plane case the lift coefficient for the N. A. C. A. 4412 
airfoil in two-dimensional potential flow is obtained as 
(cf. reference 14) 2u-(1.114) sin (a+ constant) or the 
slope of the lift curve in the monoplane case is about 
13 percent greater than that in the biplane example 
treated). Putting G L =0, 0.5, 1.0, and 1.5, respec- 
tively, the angles a e are determined as a e =— 34°51', 
— 30°8', — 25°24', and — 20°35', respectively (i. e., 
by equation (5.22) ao=-4°48', 0°5', 4°33', 9°22', 
respectively). Columns 12-15 of the table give 

the values of and for the foregoing four values 

of Oo, and are determined by equations (5.9) and 
(5.13). The velocities at each profile surface V\ and 
v t are given in terms of the stream velocity V by the 

formulas (5.23) and (5.24), and ^ = 

^2^ 2 ' Tiie P ressure ratios ^ and 2* shown in figure 

12 are given by equation (5.25). 

The Smithsonian Tables of Hyperbolic Functions 
were found useful in the numerical work. 
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TABLE I. BIPLANE ARRANGEMENT N. A. C. A. 4412 
(SEE FIG. 1Q) 
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